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We present a new theoretical method to study a trapped 
gas of bosonic two-level atoms interacting with a single mode 
of a microwave cavity. This interaction is described by a tri- 
linear Hamiltonian which is formally completely equivalent 
to the one describing parametric down-conversion in quan- 
tum optics. A system of differential equations describing the 
evolution, including the long-time behaviour, of not only the 
mean value but also the variance of the number of excited 
atoms is derived and solved analytically. For different initial 
states the mean number of excited atoms exhibits periodically 
reappearing dips, with an accompanying peak in the variance, 
or fractional collapses and revivals. Closed expressions for the 
period and the revival time are obtained. 



I. THE PHYSICAL SYSTEMS 

Since the experimental realization of weakly interact- 
ing atomic Bose-Einstein condensates confined in mag- 
netic traps 1^ much work was done to study the interac- 
tion of such a system with light. Examples include the 
scattering of light by a condensate ||^ , its refraction index 
nonlinear atom optics |^,||; the optical detection of 
the relative phase between two condensates 0, and the 
spontaneous emission in the presence of two condensates 



In this paper we are concerned with the question of how 
a condensate of two-level atoms interacts with a single 
resonant quantized cavity mode of the electromagnetic 
field. The system can be realized by placing the trap 
containing the condensate in a microwave cavity. We 
consider the limit of vanishing temperature so that es- 
sentially all atoms are in the condensate. The two states 
\g) and |e) are both members of the internal ground-state 
manifold of the condensed atoms but have a different hy- 
perfinc magnetic quantum number mp. For both states 
TO p should be chosen in such a way that the atoms remain 
trapped after the transition from \g) to |e). The internal 
energies corresponding to these states are denoted by Eg 
and Ee, respectively. Although for the existing conden- 
sates the trapping potential is different for states having 
a different magnetic quantum number we assume for sim- 
plicity that the trapping potential is the same for \g) and 



|e) . This situation may be realized by using non-magnetic 
traps or by exposing the atoms to a highly detuned laser 
beam. If the shape of the beam is appropriately cho- 
sen the induced dipole-potential can compensate for the 
difference between the trap potentials. 

Under these conditions the Hamiltonian in its second 
quantized form is given by 



i=e,g 

+hLOca'^a + Hint + Hn,i 
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where '^i, i — e,g is the field operator for atoms in the 
state \g) and |e), respectively, a is the annihilation oper- 
ator for photons in the microwave cavity. The center-of- 
mass Hamiltonian 
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is the same for both internal states. M denotes the 
atomic mass, and lu± and lOz are the trap frequencies 
in the x — y plane and in the z-direction, respectively. 
The nonlinear Hamiltonian 



Hn.i. — 
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5yVE'l(x)vI/t(f)vE-,.(f)vE-,(f) . (3) 



describes the interaction between the atoms. The quan- 
tities gij are connected to the scattering lengths for the 
scattering between atoms in the internal states \i) and 
\j). We will make the Hartree approximation and con- 
sider the case that all atoms are in the same state ipo{x) 
regarding the center-of-mass motion. This state fulfills 
the stationary nonlinear Schrodinger equation (see, e.g., 
Ref. with energy eigenvalue fj, and is the same for 
both internal states. 

We model the interaction of the atoms with the mi- 
crowave photons by the usual magnetic dipole coupling 



Hi. 



[aB + a^B*] ■ / d^x{mge^l{x)^g{x) + H.c.} 
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where rrige is the magnetic dipole moment of the atoms 
and B the magnetic field of the cavity mode at the po- 
sition of the trap. As we consider the case that the 
trap size is much smaller than the wavelength of the 
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microwave photons the interaction does not affect the 
center-of-mass motion of the condensed atoms. This jus- 
tifies the assumption that the atoms will always remain 
in the spatial mode ipo{x). 

Under these assumptions we can neglect all spatial 
modes beside fo{x) and can replace the atomic field op- 
erators vE'i(25') by ipo{x)bi{x), i = e,g, where the operators 
bi := J (P xip'^{x)'^ i{x) are the annihilation operators for 
atoms in the internal state Substituting this in the 
Hamiltonian (]l|) the latter can be simplified to 

H ^{^ + Ee)blbe + (m + Eg)blbg + hujca^a 

-[aB + a^B*] ■ {blb.rrig, + fogb^m^J . (5) 

In the interaction picture and after the rotating wave 
approximation the Hamiltonian reduces in resonance 

{huJc ^ Ee- Eg) to 

H = -hniblb^a'' + bgbla} . (6) 

The Rabi frequency f2 is given by Irhg,, ■ B*\/Ti. 

The Hamiltonian (|^) is formally equivalent to that used 
to describe parametric down conversion in nonlinear op- 
tics (see, e.g., Ref. |^). In this system 6e destroys a 
photon of the pumping beam and and 5j create a 
signal and an idler photon, respectively. The interac- 
tion between the different modes is caused by a nonlinear 
medium and describes the process that a pumping pho- 
ton is converted into two other photons under conserva- 
tion of energy. Theoretically the pumping beam is often 
considered as a classical field. To be short we will focus 
in the following on the interpretation of the Hamiltonian 
as an interaction between two-level atoms in a cavity. 
Though the physical interpretation is different our results 
can immediately be applied to parametric down conver- 
sion if one considers Ne — b\be as the number operator 
for pumping photons. 

A trilinear Hamiltonian has been examined numeri- 
cally by Walls and Barakat . Kumar and Mehta ||l|] 
derived analytical expressions for the time evolution un- 
der the condition that the number of one of the three 
particle species involved in Eq. (^) remains strongly pop- 
ulated. That the degenerate parametric oscillator may 
show collapse and revival has been predicted numerically 
by Jyotsna and Agarwal These results are similar 

to the numerical results of Drobny and Jex [T^ ] . An ex- 
act solution for a fixed number of atoms was given in a 
non-closed form by Tavis and Cummings [l4| . 

II. COUPLED DIFFERENTIAL EQUATIONS FOR 
THE MEAN OCCUPATION NUMBER AND ITS 
VARIANCE 

Our aim is to present a new analytical method for solv- 
ing the underlying operator equations. In particular this 
will then enable us to describe for the system given above 
the long-time behaviour of the mean occupation Ne{T) 



of the excited atomic mode for the case that initially all 
atoms are excited and there are zero photons in the cav- 
ity and to confront it with the time dependence of its 
variance. To do so we want to give closed-form approx- 
imate analytical solutions for Ne{T) for different initial 
states of the condensate like atomic number state, coher- 
ent state, or mixture. This will allow us to demonstrate 
that -/Ve(r) shows a periodic behaviour or the existence 
of collapse and revival, respectively. 

The Hamiltonian (0) leads to two conserved quanti- 
ties: the total number of atoms Sa '■= b\bf, + b^gbg and, 
because of the rotating-wave approximation, the num- 
ber of excitations Se '■= b\bf. -\- a^a. Because any two 
of the number operators Ne,Ng, Na can be expressed by 
the third and the two conserved quantities, it suffices to 
solve the Heisenberg equation for only one of the three 
number operators. We will consider the number operator 
Ne of excited atoms. 

Introducing the new time variable t :— fit and consid- 
ering the Heisenberg equation for as well as for 

iVe = -iblbea'' + ibgbla , (7) 

where the derivative is taken with respect to r, one easily 
finds (see, e.g., Ref. Q) 

Ne = 6iVf - 2ANe + 2B . (8) 

For notational convenience we have introduced the op- 
erators A := 2Se + 2Sa + 1 and B := SaSe- Note 
that while Sa and Se commute with every operator oc- 
curring in the problem, N^. does not commute with N^. 
This fact seems to prevent the integration of the second 
order differential equation to an equation of first order as 
it is usually done with differential equations for ordinary 
functions. A common way to circumvent this difficulty is 
the vanishing variance approximation, where a differen- 
tial equation for N^ is derived under the assumption that 
the variance Ag defined by Ag := {{N^ — Ne)'^) is neg- 
ligible. Here and in the following the mean value of an 
operator O will be represented by a bar, O :— ('0|O|'0)- 
The complete state of the system is denoted by lijj). For 
the system under consideration this approximation was 
studied by Kumar and Mehta starting from Eq. (H). 
They showed that it is valid for short times only and 
can be applied only to a very limited class of physical 
problems. 

Our objective is to improve this scheme: We want 
to go a step further and derive differential equations 
for the mean value Ng as well as for the variance Ag 
and study the respective long time behaviour for dif- 
ferent states of the condensate. Because of 9^(iVg) — 
{NeNe + NeNe + 2iV|) , thc discussiou of Ae (t) needs the 
knowledge of N^ and therefore the exact integration of 
the operator equation (H). This can indeed be done by 
observing that the commutator between N^, and A^e is 
given by [Ae, Ae] = H/(inn) so that N^N^ + NeN^ = 
{2/?,){dNl/dT) + (l/3)Ae follows, what can be used to 
integrate Eq. (||) to 
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(iV^)2 = 4Ar3 _ 2ANl + 2Ne{2B + 1) + C 



(9) 



Here C is a constant of motion which after some alge- 
bra can be written as C = —H^/{h^i}'^) + 2SaSe- It is 
therefore related to the conservation of energy and parti- 
cle numbers. Note that no approximation has been made 
up to this point. 

Although Eq. (H) cannot be exactly integrated because 
C does not commute with N^^ we are now able to treat 
together with Ae if we make the two following assump- 
tions of a vanishing asymmetry approximation: (i) the 
conserved particle numbers are always uncorrelated to 
the number of excited atoms so that (SiNji) — (5i)(A^g) 
holds for i — A, E and any integer I. This equality is ex- 
actly fulfilled for all times if the initial state is an eigen- 
state of Si, e.g., a number state, (ii) we can always ne- 
glect {{Ne — Ne)^) so that (A^e) be approximated by 
+ SA^eAg. Accordingly, in contrast to the vanishing 
variance approximation, we only assume the vanishing 
of the higher moment {{N^ — N^)^) in the hierarchy of 
moments {{Ne — A^e)™), m = 2, 3, ... of the probability 
distribution. The case m = 3 describes its degree of 
asymmetry. 

Based on this we derive from Eqs. (||) and (||) the 
system of ordinary differential equations 

dlNe ^ 6{Al + N^) - 2ANe + 2B , (10) 
d^Al = 20N^ + eOA^iVe - 8A{Al + Nj) 

+A{l + ZB)Ne + 2C -dlNl , (11) 

Eq. (|l^) may be considered as an algebraic equation for 
Ag which can be inserted into Eq. (11) to decouple the 
system. This results finally in 



dtN, 



AQAdlNe + eONed^Ne - 240N^ + 120 AN^ 
'Ne{24:-4:8B-16A^) + 12C+16AB (12) 



Ae can then be obtained from Eq. (p^. 



III. ANALYTICAL SOLUTION FOR A NUMBER 
STATE 

We are now able to treat the envisaged physical situ- 
ation: all atoms are initially for t — excited and no 
photons are present. It is desirable to obtain the devel- 
opment in time of (t) for the case of the atomic state 
being a number state with n atoms, the corresponding 
complete state of the system being denoted by \ipn)- We 
will demonstrate below that Ne{T) for other physical situ- 
ations can be reduced to this. The operator equations (||) 
and (H) show that the initial conditions for Ne{T) which 
are to be fulfilled for a number state are: iVe(O) — n, 
drNe{0) = d^Ne{0) = 0, and d^Ne{0) = -2n. 

We have in fact been able to find an exact solution of 
Eq. ([l^ ) fulfilling the correct number state initial condi- 
tions for Ne{0),drNe{0),d^Ne{0), and, to leading order 
in n, for d^Ne{0). It is given by 



7Ve™""(r) = Afe"™(0) - -muj^cn\ujT - K{m)\m) , (13) 

where the parameters m and ui are solutions of the alge- 
braic equations 

tj^m(m - l){5[6iVe"""X0) -A]- 2u?{2m - 1)} 
-120(iV^™"X0))3 60A(7V^™(0))2 ^ 8AB + 6C 
-H4Afe""™(0){3 - 2A^ - 6B} 

uj^iA + 19m(TO-l)] + 10w2(2m - 1)[A - 6Ae"™(0)] 
+ 180(iVe"""X0))2 - 60A7V™'"(0) -G + iA^ + UB . 

The function cn(z|m) is one of the Jacobian elliptic func- 
tions, and K{ra) is the complete elliptic integral |jl^. If 
m is close to one cn(z|TO) becomes almost identical to the 
expression 1/ cosh(z) except that cn^(z|m) is periodic in 
z with period 2K{m) w ln[16/(l — m)]. For a number 
state \ipn) we have A = An+1, B = , and C = 2n^ — n. 
The corresponding solution ( p3| ) is, to the first two lead- 
ing orders in n ^ 1, characterized by a; « ^/n + 2 and 
TO « l~2/n. A second exact solution of Eq. (12) is struc- 
turally similar to that of Eq. ( [T^ ) if cn(z|TO) is replaced by 
the elliptic function dn(2;|TO). In this case the parameters 
are determined by different algebraic equations which we 
will omit in this paper. 
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FIG. 1. Mean value TV""™ (solid curve) and variance 
(dashed curve) of the number of excited atoms as a function of 
time for 100 initially excited atoms in an atomic number state. 
The dot-dashed curve is obtained in the vanishing variance 
approximation, fl is the generalized Rabi frequency. 

Our improved mean value solution Afe™(T) of Eq. (|l|) 
is plotted in Fig. 1 together with the ordinary mean value 
solution obtained in the vanishing variance approxima- 
tion (Ae =0). As one can see, the fluctuations become 
very large when the dip in the number of excited atoms 
occurs. As compared with the vanishing-variance ap- 
proximation the dip is reduced by a factor of 1/2. This 
is in qualitative agreement with the results of Ref. |l^ ] 
where a similar physical system was examined under the 
assumption that the photons escape quickly so that the 
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number of excitations is not conserved. Note that the fre- 
quency of the occurrence of the dip is modified, too. The 



period Tp is determined by ^/n 
leading to Tp = \n{8n) / [^/r^- 
increasing number of atoms. 



2nTp = 2K{1- 2 In) 
2r2]. It diminishes with 



IV. COLLAPSE AND REVIVAL 

Based on the solution (|l^ we are now able to cal- 
culate the mean number of excited atoms if initially no 
photons and only excited atoms are present which form 
a coherent state with mean number h of atoms. We ex- 
pand this state in the number states |'(/'/). Since each 
of the states •\/iVe(0) exp[— is an eigenstate of 
Sa with eigenvalue / they are orthogonal for different 
This reduces the mean value of N^^ to the expression 



(=0 



(14) 



Because the r.h.s. is only a function of the mean values 
iV""™(T) of the number state case, it may directly be 
evaluated with the help of Eq. ( [l3| ) . The physical reason 
for this simplification is that the interaction conserves 
the total number of atoms. Thus, no interferences be- 
tween number states of different total atom number can 
occur. Note that the expression is identical to the 
mean value for a statistical mixture of number states with 
Poissonian statistics. 




FIG. 2. For a coherent state of 100 initially excited atoms 
or a Poissonian mixture of number states the time evolution 
of the mean number of excited atoms exhibits collapses and 
revivals. 

Whereas 7V""™(t) for a number state shows a periodic 
behaviour with period Tp it is an important consequence 
of Eq. ( p^ that a coherent atomic state as initial state 
leads to collapse and revival in the number N'^^ij) of ex- 
cited atoms (see Fig. 2). This is also the case for a statis- 
tical mixture. Leaving open the question if a condensate 



can at all be described by a coherent state (compare, e.g., 
Ref. ||r^]) an experimental verification of collapse and re- 
vival cannot be used to discriminate between coherent 
state and mixture. 

We turn to the details of the time development of 
N'^°^(t). The time Tr when the first (large) revival oc- 
curs can be determined by the criterion given in Ref. p8[ : 
the expression (V'i|iVe(''')|V'i) must simultaneously have 
a maximum for the two neighboring values / = n and 
I = n -|- 1. This guarantees that many terms in the 
sum in Eq. ( p^ ) become simultaneously large. Since 
cn(z|m) is periodic this condition becomes a condition 
on the argument of the Jacobian elliptic function, namely 
VnnTR-K{l-2/n) = 2rK{l~2/n) and y/n + T^TR- 
K{1 - 2/{n + 1)) = 2(r -t- l)K{l - 2/{n -h 1)) with an 
integer r. This can be solved for r and Tr and results to 
leading order in n in 



Th 



2^/fl \n{%n) 
9. ln(8n) - : 



(15) 



According to its derivation the revival time Tr is inde- 
pendent of the probability weights of the mixture or state 
when it is expanded in | ?/';). It is therefore the same for 
all mixtures with the same initial mean number of excited 
atoms. 

In addition to the large revival at time Tr there are 
many smaller revivals at earlier times which arise at frac- 
tions of Tr (see Fig. 3). These fractional revivals have 
been predicted numerically for the degenerate paramet- 
ric oscillator in Refs. W2 



The fractional revivals arise 
when only a part of the sum terms in Eq. ( p^ have si- 
multaneously a maximum. Their revival times can be 
calculated by assuming that the terms for I = n and 
I = n-\- r, r = 2, 3, . . . are simultaneously maximal. This 
leads to TR{r) — TR{\)/r. Since the revivals reappear pe- 
riodically one generally finds them at fractions of small 
integers of the time Tr. 




FIG. 3. Between the collapse and the main revival frac- 
tional revivals do appear. 
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An estimation of the almost constant value iVconst 
taken by iV™^(r) between the collapse and the first re- 
vival can be derived under the assumption that the argu- 
ments of cn(z|TO) of Eq. (|l^ ) that appear in the sum in 
Eq. are statistically distributed over the whole pe- 
riod of cn(z|m) (i.e., z e [0, 2/-ir(m)]). A^const then is sim- 
ply the average of Eq. (O) over one period Tp. The corre- 

sponding integration {2K)~^ cii^{z\m)dz can be per- 
formed if cn^(z|m w 1) is approximated by 1/ cosh^(z) 
and results to leading order in n in 

iV_.«..{l-^}. (16) 

The two analytical expressions for Tr and iVconst ai'e 
in good agreement with the numerical evaluation of 
Eq. (0). 

In conclusion we have presented analytical results for 
the collapse and revival of the mean occupation num- 
ber and the variance of excited atoms in a microwave 
cavity. The phenomenon of a revival in the context of 
Bose-Einstein condensation was also examined for differ- 
ent physical systems. Zhang and Walls |]l9| discovered it 
numerically for atoms passing through a standing light 
wave. Wright et al. showed that it occurs if the Bose 
condensate is described by a superposition of states with 
different total number of atoms. Our results differ from 
these approaches in that a third mode (the microwave 
photons) is fully incorporated and that the new theoret- 
ical method applied here allows to derive closed expres- 
sions for the interesting physical quantities. Acknowl- 
edgement: We thank Markus Holzmann for discussions 
and the Optik-Zentrum Konstanz for financial support. 
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